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Abstract. I review recent progress in developing a complete dynamical model for the evolution 
of the Little Bang fireballs created in relativistic heavy-ion collisions, and using the model to 
Ft extract the transport properties and initial density fluctuations of the liquid quark-gluon plasma 

' . ' state of matter of which makes up these Little Bangs during the first half of their lives. 

<N 

R -1 1. The Big Bang and the Little Bangs 

Ultra-relativistic heavy-ion collisions at RHIC and the LHC produce fireballs made of 

t-h extraordinarily hot matter, at initial energy densities (at the time when the matter reaches 

approximate local thermal equilibrium) that exceed the energy density of atomic nuclei in their 

,—1 ground states by two to three orders of magnitude. Due to enormous pressure gradients between 

<£ the fireball center and the surrounding vacuum, these fireballs undergo explosive collective 

expansion, cooling down rapidly through several different states of matter, finally fragmenting 

into thousands of free-streaming hadrons whose energy and momentum distributions can be 

detected in the detectors set up around the collider rings. The evolution history of these 

"Little Bangs" has much similarity with the Big Bang that created our Universe (Fig. [T]): Both 

undergo Hubble-like expansionjj feature a hierarchy of decoupling processes that are driven by 

the expansion dynamics (and not by the finite geometric size of the Little Bang), with chemical 
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decoupling of the finally observed particle abundances 



initial-state quantum fluctuationsjj (shown in Figures 
experimentally observed final state. 
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1 The relative velocity between two matter elements increases roughly linearly with their relative distance 
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2 In the Big Bang the chemical composition is frozen during the process of primordial nucleosynthesis at an age 
of about 3 minutes, in the Little Bang this happens during hadrosynthesis at the quark-hadron transition at an 
age of about 20-30 ioctoseconds [5]. 

3 In the Big Bang, the formation of neutral atoms by charge recombination at an age of about 380,000 yr makes the 
universe transparent to light and freezes the thermal Bose- Einstein energy distribution of the Cosmic Microwave 
Background radiation. When the Little Bang reaches an age of about 40-50 ioctoseconds, the final-stage hadron 
gas becomes so dilute that strong interactions between the hadrons cease and their energies and momenta are 
"frozen out" . 

4 The initial wave function of the universe, stretched by cosmic inflation, seeds density fluctuations in the Big 
Bang. In the Little Bang, the initial quark and gluon wave functions inside the nucleons within the colliding 
nuclei control its initial density distribution. 

5 In the Big Bang the initial density fluctuations evolve under the action of gravity, described by Einstein's general 
theory of relativity, into the observed Cosmic Microwave temperature fluctuation spectrum shown in Figs. [2] (left) 
and[3](top), and ultimately into today's distributions of stars, galaxies, and galaxy clusters and superclusters (see 
top panel of Fig.fTh. In the Little Bang, they evolve through viscous hydrodynamics into the measured anisotropic 
collective flow patterns and their event-by-event fluctuations that are the subject of this overview. 
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Figure 1. Artist's conception of the evolution of the Big Bang (top 
(bottom - credit: Paul Sorensen and Chun Shen). 



credit: NASA) and the Little Bang 



Of course, the Big and Little Bangs are quite different in other aspects: Their expansion rates 
differ by about 18 orders of magnitude; the Little Bang's expansion is 3-dimensional and driven 
by pressure gradients, not 4-dimensional and controlled by gravity; Little Bangs evolve on time 
scales of ioctoseconds, not billions of years; distances are measured in femtometers rather than 
light years. Most importantly, the Little Bang Standard Model is still under construction. This 
overview discusses recent progress of the edifice. 



2. Eccentricity fluctuations, anisotropic flows, and flow fluctuations 

We can observe only one Big Bang (the one that produced our universe), but at the Relativistic 
Heavy-Ion Collider (RHIC) and Large Hadron Collider (LHC) we have experimentally created 
and studied billions of Little Bangs. Each Little Bang is different: Highly successful 
phenomenology based on hydrodynamic evolution models Uj has taught us that the initially 
very dense quark-gluon matter created in heavy-ion collisions reaches approximate local thermal 
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Figure 2. Temperature fluctuation spectrum of the Big Bang at age 380,000 yr (left), as measured through the 
Cosmic Microwave Background radiation by WMAP, and of four typical Little Bangs created in central Pb+Pb 
collisions at the LHC (right), at age 0.2fm/c, as calculated from the IP-Glasma model [T]. Figure taken from [2]. 



equilibrium on a very short time scale [9j of order 1 fm/c (3 ioctoseconds) , after which it evolves 
according to the macroscopic laws of relativistic viscous fluid dynamics. Pressure gradients in 
the fluctuating density profile (see Fig.Huor a specific example) are the hydrodynamic forces that 
accelerate the fluid and cause it to expand and dilute. Spatial anisotropics and inhomogeneities 
in the initial density profile transverse to the beam direction lead to corresponding anisotropics 
in the final transverse expansion flow velocity profile that are imprinted on the momenta of 
the experimentally observed particles emitted from the collision. Each Little Bang features 
its own final flow velocity profile which (within hydrodynamics) is a deterministic classical 
response to the initial conditions that fluctuate from collision to collision due to quantum 
fluctuations in the initial nuclear wave function. We now know that the quark-gluon liquid 
that makes up the matter of the Little Bang during the first half of its life has very small 
viscosity, behaving like an almost ideal fluid [8]. This is a fantastic gift of Nature since it 
allows us to study experimentally the spectrum of initial-state quantum fluctuations through 
the final-state anisotropic flow fluctuations. Had the quark-gluon plasma (QGP) turned out to 
be highly viscous, all initial-state fluctuations would have been wiped out by dissipation before 
final decoupling of the emitted particles, thereby closing this observation window on the initial 
state of the Little Bang and on the quantum nature of the initial energy deposition process. 

The Little Bang pressure p is related to its energy density e by its equation of state (EOS), 
and the anisotropy of the initial pressure gradients can be characterized by a series of harmonic 
eccentricity coefficients 
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where e(r, p) is the initial energy density distribution in the plane transverse to the beam 
direction. e n characterizes the magnitude of the nth harmonic deformation coefficient, and the 
angle ^f n gives the direction of the corresponding deformation component in the lab frame. Both 
the magnitudes e n and their orientations *$> n fluctuate from event to event. The mean values 
{e n ) averaged over many events from three popular initial energy deposition models (see [1] for 
a description and original references) are shown as a function of harmonic index n in the left 
bottom panel of Fig. |3| They represent the primordial temperature fluctuation spectrum of the 
Little Bang. We see that each collision centrality generates a different class of Little Bangs, 
with its own (e n ) power spectrum, and that the three models shown give quite different results 
for these power spectra. If we can somehow measure the centrality dependence of the initial 
(e n ) power spectrum, we have a powerful constraint on the initial nuclear quark and gluon wave 
functions on our hands. This is reminiscent of the crucial role the CMB temperature power 
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Figure 3. Top: Temperature power spectrum of the Big Bang at age 380,000 yr, as measured by the Planck 
satellite 3 . Bottom left: Primordial eccentricity power spectrum of the Little Bangs created in 2.76 A TeV 
Pb+Pb collisions of different centralities, at t — 0, from the energy density distributions of three different initial- 
state models (IP-Glasma, MC-Glauber, MC-KLN) 4,. Bottom right: The final Little Bang flow power spectrum 
for ultracentral (0-0.2% centrality) Pb+Pb collisions at y/s = 2.76 A TeV at the LHC, measured by the CMS 
Collaboration (Wei Li, Quark Matter 2012) and calculated with viscous hydrodynamics VISH2+1 [6] using MC- 
Glauber and MC-KLN initial conditions with the indicated specific shear viscosities. Figure taken from [7]. 



spectrum (shown in the top panel of Fig. ufy has played in nailing down the parameters of the 
Standard Model of Big Bang Cosmology. 

The hydrodynamic response to the initial e n -spectrum of the Little Bang is reflected in the 
observed transverse momentum distributions dNi/ {dypTdpTdcj)) of the various emitted hadronic 
species i, and in correlations among them. The transverse momentum spectra can again be 
characterized by a set of complex harmonic coefficients 
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Figure 4. Typical transverse energy density profiles e(x, y) from the IP-Glasma model [T] for a semiperipheral 
(6 = 8 fm) Au+Au collision at y/s = 200 A GeV, at times r = 0.01, 0.2, and 5.2 fm/c. From r = 0.01 fm/c to 0.2 fm/c 
the fireball evolves out of equilibrium according to the Glasma model [10J : at r = 0.2 fm/c the energy momentum 
tensor from the IP-Glasma evolution is Landau-matched to ideal fluid form (for technical reasons [IT] the viscous 
pressure components are set to zero at the matching time) and henceforth evolved with viscous Israel-Stewart 
fluid dynamics, assuming r//s — 0.12 for the specific shear viscosity. The pre-equilibrium Glasma evolution is seen 
to somewhat wash out the large initial energy density fluctuations. The subsequent dissipative hydrodynamic 
evolution further smoothes these fluctuations. The asymmetric pressure gradients due to the prominent dipole 
asymmetry in the initial state of this particular event (visible as a left-right asymmetry of the density profile in 
the left panel) is seen to generate a dipole ("directed flow") component in the hydrodynamic flow pattern that 
pushes matter towards the right during the later evolution stages. 



Here <f> is the azimuthal angle around the beam direction of the particle's transverse momentum 
PT) and the curly brackets denote the average over particles from a single collision. Each particle 
species i has its own set of anisotropic flow coefficients V n . Eq. ([2]) defines the flow coefficients 
and associated flow angles for the entire event, whereas Eq. (13J) is the analogous definition for 
the subset of particles in the event with a given magnitude of the transverse momentum pr- I 
suppress the dependence of both types of flow coefficients on the rapidity y. v n are known as the 
"integrated" anisotropic flows, v n (pr) are called "differential" flows. By definition, both v n and 
v u(pt) are positive definite. Hydrodynamic simulations show that in general the flow angles Vl/ n 
depend on pj, and that, as a function of pr, ^ n {pr) wanders around the "average angle" VP n 
that characterizes the integrated flow v n of the entire event [12] . 

For each collision system, centrality class, and collision energy, the fluctuating initial state of 
the Little Bang is characterized by a distinct probability distribution P(e n ,$ ri ), one particular 
moment of which is the (e n ) power spectrum shown in Fig. K3l For different n, the angles 
3> n are more or less correlated with the direction of the impact parameter b, but this direction 
cannot be measured experimentally. The viscous hydrodynamic evolution relates the probability 
distribution for the initial complex eccentricity coefficients deterministically to probability 
distributions P(V n ) and P(V n (pr)) for the integrated and differential final complex harmonic 
flow coefficients; these distributions are particle species specific. As will be discussed in the next 
section, the relation between the final complex flow coefficients V n and the initial eccentricity 
coefficients e n e m * n (and thus between the corresponding flow and eccentricity probability 
distributions characterizing each class of collisions) depends on the viscosity of the Little Bang 
matter. One goal of the relativistic heavy-ion program is to both constrain the QGP viscosity and 
identify the correct theory for the initial-state quantum fluctuations by performing a complete 
experimental reconstruction of the final multi-dimensional distributions P(V n ) and P(V n (pr)). 

Due to limited statistics arising from the finite number of particles emitted by each Little 
Bang, neither the magnitudes v n nor the flow angles ^ n can be accurately determined for a single 
event. Experimental flow measures therefore involve angle correlations between two or more 
particles; for example, ({e m</il }p T1 {e- m<fc }p T2 ) = (v n (p T i)v n (p T 2) cos[n(^ n (p T1 )-^ n (p T2 ))]} 
where the first particle from the event has transverse momentum pjl = (pti > 4>i ) and the second 
particle has PT2 = (PT2, §2 ), the average {. . .} is over all such particles in the event and (...) 
indicates the average of the result over many Little Bangs of the selected class. One sees that 
the event-by-event fluctuations and transverse momentum dependences of both the anisotropic 
flow magnitudes v n and their associated flow angles *& n affect these experimental observables. 



Different such observables correspond to different correlation functions between the v n 's and 
^>nS all of which can, for a given initial energy deposition model with probability distribution 
P(e n ,$n), be computed from the hydrodynamically predicted probability distributions P(V n ) 
and P(V n (pr)). One can define experimental observables that separate the fluctuations of the 
anisotropic flow magnitudes v n (pr) from those in the flow angles ^ n (pr) [12] • These should be 
powerful discriminators between different P(V n ) and P(V n (pT)) distributions, and thus between 
different initial-state fluctuation models. Experimental studies of these v n and ^f n fluctuations 
and their pp-dependences have just gotten under way; the results will be interesting and should 
significantly advance the construction of the Little Bang Standard Model. 
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Figure 5. Three-plane correlations between the initial participant planes (left 6-panel) and the final-state 
flow planes (right 6-panel), for 2.76 ^4TeV Pb+Pb collisions as functions of collision centrality, indicated by the 
number of participant nucleons -/V part on the horizontal axes. The solid (dashed) lines are for MC-Glauber (MC- 
KLN) initial conditions evolved event-by-event with viscous hydrodynamics using jj/s = 0.08 (0.2) for the specific 
shear viscosity [14]. Filled circles show the experimental values measured by the ATLAS Collaboration |16j . 

The angles <3? n associated with the initial deformation parameters e n are not only correlated 
with the (unmeasurable) direction of the impact parameter b, but also with each other, in ways 
that are predicted by (and thus depend on) the initial energy deposition model [HI Q3]. These so- 
called participant-plane correlations in the initial state are translated by hydrodynamic response 
into final-state flow angle correlations. These flow angle correlations were measured in Pb+Pb 
collisions at the LHC by the ALICE [15] and ATLAS [E] Collaborations. All measured two- and 
three-plane flow angle correlations [16] are qualitatively reproduced by viscous hydrodynamic 
calculations, both in their magnitudes and centrality dependences [TJ] (right panels in Fig. M. 
They differ, however, qualitatively from the initial participant-plane correlations [H], often even 
in sign (left panels in Fig. foj . Hydrodynamic evolution is nonlinear and leads to mode-coupling 
between different harmonics [17J. For example, elliptic and triangular deformations £2 and £3 
provide a nonlinear contribution to pentangular flow V5 which, for large impact parameters 
where the elliptic deformation £2 is big, overwhelms the linear response to £5 and completely 
decorrelates the pentangular flow plane ^5 from the angle $5 of the initial pentangular density 
deformation [T7]. Only by accounting for mode-coupling, either directly by following the non- 
linear hydrodynamic evolution [14] or through a non-linear response analysis that keeps at least 
second-order terms |18j . can the experimental data be reproduced. These flow-plane correlations 
thus represent an experimentum crucis in support of the hydrodynamic paradigm for the Little 
Bang; dynamical models without a large degree of local thermalization and hydrodynamic 
collective flow will not be able to describe the dynamical change of character between the initial- 
state participant-plane and final-state flow-plane correlations n " " 



The strength of the flow-plane correlations depends on the QGP shear viscosity [14 . A careful quantitative 



3. QGP shear viscosity from anisotropic flow measurements 

The efficiency of the fluid to convert spatial anisotropies in the pressure gradients into anisotropic 
flows is degraded by shear viscosity. The "conversion efficiency" v n /e n is therefore a measure 
for the specific shear viscosity rj/s of the expanding fluid. This is shown in Fig. [6J Higher 
harmonics, reflecting variations on smaller spatial scales, are suppressed more strongly by shear 
viscosity than lower flow harmonics. 
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Figure 6. Left: The eccentricity-scaled integrated elliptic flow of all charged hadrons, u| : (v I ' s ) / ' £ ■> as a 
function of total charged hadron multiplicity density per unit overlap area, (l/S)(dN c h/dy). The experimental 
data points show two measures for the elliptic flow ((i>2) [S] and «2{2} [20]) from 200ylGeV Au-Au collisions at 
RHIC, measured by the STAR Collaboration. Both panels use the same sets of data, but use different average 
initial eccentricities (e) and overlap areas (S) (obtained from the MC-Glauber and MC-KLN models) to normalize 
the vertical and horizontal axes. The theoretical curves were computed with the VISHNU model [21], for different 
(temperature-independent) choices of the specific QGP shear viscosity (jj/s)qgp- Right: Viscous suppression of 
V n for T]/s — 0.08 (squares) and 0.16 (circles), as function of harmonic index n, for Au+Au collisions at 20-30% 
centrality |22| . 



The azimuthally symmetric part of the collective transverse flow, called radial flow, boosts 
particles from low to high transverse momentum in proportion to their masses; it is thus 
responsible for the distribution of the hydrodynamically generated momentum anisotropies over 
the various particle species and in px [2H]. To correctly describe the differential anisotropic 
flows v n (pr), for all charged hadrons or for specific identified hadron species, thus requires 
that the total hydrodynamic momentum anisotropy (which "measures" rj/s), the radial flow 
(which is also sensitive to bulk viscosity that suppresses radial flow |26j). and the chemical 
composition of the system at final decoupling (which reflects the kinetics of chemical freeze-out) 
are all correctly described by the model. By analyzing instead the total flow anisotropies i^ h 
of all charged hadrons, integrated over pr, one can strongly reduce the model sensitivity to 
bulk viscosity and final chemical composition and obtain a much more robust estimate of rj/s 
|21| I25j . The left panel in Fig. p\ shows such an extraction of rj/s from elliptic flow data. The 
VISHNU model used for this extraction couples the viscous hydrodynamic evolution of the QGP 
phase to a microscopic kinetic evolution of the hadronic phase after hadronization, thereby 
eliminating many uncertainties in earlier simplified calculations that described the hadronic 
phase macroscopically, too. The different theoretical lines show that, for a given initial energy 
deposition model, vf'/£2 is only sensitive to the specific QGP shear viscosity. Once (?7/s)qgp 
has been adjusted to the measured total charged hadron elliptic flow as in Fig. p\ (resulting in 
(rj/s)QGP — 0.2 for MC-KLN and (?7/s)qgp ~0.08 for MC-Glauber initial conditions), one finds 
that the model also correctly describes the pr- s P ec tra and differential elliptic flows V2(pt) of all 
charged hadrons together, as well as for specific identified hadron species (pions, kaons, protons), 
for all collision centralities [22]- Furthermore, it correctly predicted the analogous observables 
for Pb+Pb collisions at the LHC [HI E3 [2H] , in particular an increased mass splitting between 
the differential elliptic flows for light and heavy hadrons [6] [29] . 

study of these correlations by more precise experiments and more systematic theoretical analyses should be able 
to separate the effects arising from the initial-state fluctuation spectrum and from dissipative transport. 
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Figure 7. Top left: The centrality dependence of v n {2} from 2.76 ^4 TeV Pb+Pb collisions measured by 
ALICE [15] compared to viscous hydrodynamic model calculations 11 . Bottom left: Comparison of v„(pt) for 
the same collision system at 20—30% centrality from ATLAS [5H] with hydrodynamical calculations, using both 
a constant average and a temperature dependent rj/s QT]. Right: Scaled distributions of 1)2,3,4 (from top to 
bottom) from viscous hydrodynamics with IP-Glasma initial conditions [11] compared with experimental data 
from ATLAS [24] and with the scaled distributions of the corresponding initial eccentricities £2,3,4- Nonlinear 
hydrodynamic evolution causes slightly larger variances for the v n distributions compared to those of e n . The 
data in both panels are from Pb+Pb collisions at the LHC. 

Unfortunately, the two left panels in Fig. [6] also show that different initial energy deposition 
models lead to very different estimates for (j]/s)qgp, due to their different initial eccentricities 
£2 (see Fig. k3j). In [30] it was shown that this model ambiguity can be resolved by analyzing 
simultaneously elliptic and triangular flow (v2 and V3) data. This analysis eliminated the MC- 
KLN model with (t]/s)qgp = 0.2 as a viable candidate, and thereby demonstrated the power of 
a comprehensive set of anisotropic flow data to (over) constrain the dynamical evolution model. 
Unfortunately, the MC-Glauber model did not enjoy a much longer life itself: The v n data 
(n = 2, . . . ,7) from ultra-central Pb+Pb collisions shown by CMS at Quark Matter 2012 (see 
the right bottom plot in Fig. K3| and the event-by-event v n probability distributions (n = 2, . . . , 4) 
measured by ATLAS and shown at the same conference [21] (see the right panel in Fig. u\ can 
not simultaneously be described by viscous hydrodynamics for any choice of (??/s)qgp if initial 
fluctuation spectra from either of these two models (MC-KLN and MC-Glauber) are used. 

Fortunately, a new and much more successful initial-state model saved the day: the IP- 
Glasma model [311 IT], based on the Color Glass Condensate idea [10], implements gluonic field 
fluctuations inside nucleons [33|IT] as well as gluon saturation effects [10|I32], As shown in Fig. Vl\ 
these IP-Glasma initial conditions (represented by the solid lines in Fig. pi lower left panel), 
when (after a short initial pre-equilibrium stage modeled by classical Yang-Mills evolution) 
evolved with viscous fluid dynamics, reproduce the entire measured spectrum of charged hadron 



anisotropic flow coefficients v n , both integrated over and differential in p^, for all collision 
centralities, as well as the measured (23] event-by-event distribution of i>2, t>3, and V4, again 
for a range of collision centralities |11| . This is shown in Fig. u\ for Pb+Pb collisions at the 
LHC but also holds for Au+Au collisions at RHIC [UJ. The only difference is that for LHC 
energies the effective specific shear viscosity of the QGP fluid must be chosen somewhat larger 
(( ? ?/ s )lhc — 0.2) than at top RHIC energy ((??/s)rhic — 0.12). Both data sets are consistent with 
a temperature dependent specific shear viscosity (rj/s)(T) proposed in [34] (the thin lines in the 
lower left panel in Fig. FFI) that has a minimum value around l/(4-7r) = 0.08 at the pseudocritical 
temperature T cr where the QGP hadronizes and then rises to about 5 times that minimal value at 
2T cr . While this temperature-dependence of rj/s is not yet tightly constrained by the available 
analyses of RHIC and LHC data, the need for a somewhat larger effective shear viscosity at 
the higher temperatures probed at the LHC, already noted in earlier work |35t I34j. appears to 
solidify, and the early evidence from RHIC data for a very low specific shear viscosity near T cr 
[36] . perhaps as low as the KSS bound [37J of l/(47r), is strongly supported by the work in [IT] . 



4. Summary: status of and prospects for the Little Bang Standard Model 

The Little Bang Standard Model is still under construction, but its key features are showing 
through the scaffolding: 1. Every heavy-ion collision system, centrality class, and collision 
energy generates a different class of Little Bangs, each with its own characteristic initial 
density fluctuation and final anisotropic flow fluctuation spectrum. 2. The initial fluctuation 
spectrum can be computed from the Color Glass Condensate (CGC) theory using e.g. the IP- 
Glasma model; gluon field fluctuations inside the nucleons within the colliding nuclei play an 
essential role in this spectrum. 3. After a very short pre-equilibrium evolution stage, not lasting 
much longer than 1 fm/c, which is best described using classical Yang-Mills field dynamics 
corrected for quantum fluctuations, the Little Bangs are in the QGP phase which undergoes 
viscous hydrodynamic evolution until it has cooled down to the pseudo-critical temperature 
for hadronization, T cr ~160MeV. The specific QGP shear viscosity is small, around l/(4-7r) 
near T cr , rising modestly at higher temperatures; its effective value for 200^4GeV Au+Au 
collisions at RHIC is around 0.12 = 1.5/(4-7r), for 2.76 ATeV Pb+Pb collisions it is around 
0.2 = 2.5/(47r). This small shear viscosity allows flow anisotropies to build in response to the 
initial fluctuations and anisotropies in the pressure gradients, to values that are large enough to 
be experimentally observed in the final state but still visibly attenuated by shear viscous effects, 
following an ?7/s-dependent pattern where higher harmonic flow coefficients are suppressed more 
strongly than lower harmonics. The spectrum of flow anisotropies, their magnitudes, directions, 
particle species and px-dependences, and event-by-event fluctuations provide a rich menu of 
experimental observables from which the QGP shear viscosity and initial fluctuation spectra 
can be reconstructed. 4. After hadronization, the Little Bang continues to evolve as a dilute, 
highly dissipative hadron resonance gas that is highly inefficient in generating any additional 
momentum anisotropies but also doesn't erase those established earlier during the QGP phase. 
This makes it possible to extract from hadronic final state observables quantitative information 
about the earlier QGP and CGC phases. 

In this relatively simple form, the Little Bang Standard Model applies to Little Bangs created 
at top RHIC and LHC energies. At lower collision energies, our understanding of the initial 
conditions becomes less reliable and the lifetime of the hydrodynamic QGP stage shrinks, making 
a quantitative theory of the Little Bang more challenging, but not less interesting. 
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